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Abstract 

Existence of thermodynamically stable composite vortices in a two-component superconductor may form 
distinctive vortex patterns, and may lead to type- 1.5 superconductivity. Here we study the surface energy of 
the two-component superconductor and show that the sign of surface energy is determined not only by the 
Ginzburg-Landau parameters /c, (/ = 1,2) of two superconducting components, but also by a temperature 
independent parameter k^, which is defined as the ratio of the coherence lengths of two components. Since 
the negative surface energy conduces to the invasion of thermodynamically stable composite vortices into 
a superconductor, the criterions for stability of composite vortex are these three independent dimensionless 
parameters. We find that there can exist thermodynamically stable composite vortex in a type-l-i-type-2 
or type-2-i-type-2 material. We also predict that unusual vortex patterns hke those observed in MgB2 (V. 
Moshchalkov et al, Phys. Rev. Lett., 102, 117001 (2009)) can occur in some type-2-i-type-2 superconduc- 
tors. 

PACS numbers: 74.25.Ha, 74.25.0p, 74.20.De 
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It is well known that a conventional superconductor can be categorized as type-1 or type-2, 
depending on its behavior under a magnetic field. In a type-1 superconductor, such as lead or 
aluminum, superconductivity loss happens suddenly as the field surpasses a critical level. A type- 
2 superconductor, such as niobium, has two critical fields, and if the applied field is stronger than 
the lower but weaker than the higher, then the field can penetrate the material by vortices [1]. The 
criterion that determines whether a superconductor is of type-1 or type-2 is the Ginzburg-Landau 
(GL) parameter kM\. It is defined as the ratio of the penetration depth A over the coherence length 
^, K = A/^, and the critical value Kc = 1/ V2 represents the demarcation line between type-1 
(k < Kc) and type-2 (a" > Kc) superconductors. Interaction between vortices is attractive in a type-1 
superconductor Isl,^] and repulsive in a type-2 superconductor [5]. 

Novel phenomena may be found in the newly discovered two-component superconductors, such 
as MgB2 M\- In the experiment by Moshchalkov et al, unconventional stripe- and gossamerlike 
vortex patterns have been directly visualized by Bitter decorations on high quality MgBj single 
crystals OD. These observations are attributed to the coexistence of n and a superconducting 
components in MgB2. These two components are in different regimes: Kj; = A„l^n = 0.66 < 1 / V2 
(type-1) and Ko- = Ag-f^a- = 3.68 > 1/ V2 (type-2). There are evidences of existence of composite 
vortex, which is axisymmetric, and has a normal core with phase of each condensate changes In 
around the core. Interaction between these composite vortices is short-range repulsive and long- 
range attractive. This is different from that of type-1 or type-2 superconductivity, and leads to 
type- 1.5 superconductivity tAm- 

The key features of type- 1.5 superconductivity are: (i) existence of thermodynamically sta- 
ble composite vortex; (ii) interaction between composite vortices is short-range repulsive and 
long-range attractive. It has been shown in numerical simulations that the composite vortex in 
a type-l-l-type-2 superconductor is thermodynamically stable when the disparity between coher- 
ence lengths ^1 and ^2 is extremely large [8]. However we need definite criterion for stability of 
composite vortex in a two-component superconductor. 

Motivated by recent interest in exotic vortex state in multi-component superconductors, here 
we study one-dimensional superconducting-normal boundary in an abstract two-component sys- 
tem. We show that, the sign of surface energy in such system is determined not only by the 
Ginzburg-Landau parameters /c, (/ = 1, 2) of two superconducting components, but also by a tem- 
perature independent parameter k^, which is defined as the ratio of the coherence lengths of two 
components: = ^1/^2- Since the negative surface energy conduces to the invasion of thermody- 
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namically stable composite vortices into a superconductor, the criterions for stability of composite 
vortex are these three independent dimensionless parameters. The same method has been used to 
speculate the stability of vortex in conventional one-component superconductor [5]. The vortex 
discussed in present work is axisymmetric composite vortex with normal core, and phase of each 
condensate changes 2n around the core. Of course, in principle, there can exist other kinds of 
composite vortices unlike those we discussed here. But it was shown that the composite vortex 



discussed in current work is energetically preferred in weak external field [|7|,l8|,|9D. The invasion of 
composite vortex is thermodynamically favorable for certain values of Ki,K2, k^, which ensure the 
negative surface energy. We find that there can exist thermodynamically stable composite vortex 
in a type-l+type-2 or type-2+type-2 material. 

We start with a general system in which two superconducting components coexist. The GL free 
energy density of the system is 

i= 1 ' 

where /„o is the free energy density of the body in the normal state in the absence of the mag- 
netic field, Vd^PiP) = -I- bi\^if/2 (i = 1,2). 7/ is a coefficient characterizes Josephson 
coupling between two superconducting components. In the following we do not consider cou- 
pling effect and set 77 = 0. We also assume that the effective mass m* and charge e* of two 
components are equal: m* = m*, e* = e*. There are four characteristic lengths: the penetra- 
tion depth Aj and coherence length ^, for each component are given by: At = (m*c^ lAne*^'^^^)^^^, 

= h/(2m*\ai\y^^, where *F,o = {-UilbiY^^. The thermodynamic critical magnetic field of 
the individual component is Hc,{i) = Oo/(2 V5;r/l,^;), where Oq = hc/e* is the flux quantum. 
The magnetic field penetration depth and the thermodynamic critical magnetic field of the sys- 
tem © are: i = (l/Aj + l/AjT^'^, H^t = + Hl^^^Y'^. Notice that A < min(^i,i2), 
Hct > max(//rt(i), 7/^(2))- 

The problem of the thermodynamical stability of the composite vortex with normal core can 
be converted into the surface energy problem. The stability of the composite vortex depends upon 
the surface energy, or, more accurately, upon the sign of surface energy. Let us consider a plane 
interface between normal (n) and superconducting (5) phases in a two-component superconductor, 
taking the interface as the jz-plane and the x-axis into the s phase. Surface energy ans is defined as, 
under the thermodynamic critical magnetic field =HctZ, the Gibbs energy difference between 
the n, s transitional state and the fully normal state (or fully superconducting state since these must 



be equal) of the superconductor with unit cross-section: 

^■^1 2-r U l(V - + y(|VFi (H,, - V x A)H . (2) 

The integrand vanishes, both within the n phase (x -co), where = and V x A 
and within the s phase (x oo), where = *P,o and V x A =0. Now the distribu- 
tion of all quantities depends only on the coordinate x. This fact enable us to choose 
gauge potential as A = (0,Ay(x),0). Then the order parameters *F, can be taken real. We 
shall use the dimensionless quantities as following: p = x/A, if/i = ¥i/*Pio, 1/^2 = 
*F2/*I'2o. ^ = |A| IHctX, K = B = |V X A| /Z/^- Then the expression Q becomes Qr„, = 
(H^^Am II dp [ZL (Hctu)/H,,)^ [2(^,/^)V? + (H,JH,r(i^nA/AifA^i/f^ - lip] + i^f] + (A' - If] 
It can be verified that all coeflicients in the integrand can be represented as functions of three 
dimensionless temperature independent parameters: ki = A\l^\, K2 = A2/^2, = ^1/^2 llOO- And 
the surface energy ^ can be rewritten as: 
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where Ai = l/zCj+Af^/Afj, Bi = {Kl+KlKp/iKl+K^iKp, Ci = Kl/iKj+K^Kp, A2 = 1 1 kI+1 I k]k^^, B2 = 
{k\ + k\kPI[k^^{k\ + /fj/c^)], C2 = kIk^/(kI + /Cj/c|). The GL eqs of motion following from the free 
energy ([T]) are: 

A21A2 = y^V2 - <A2 + iA2' 

a" = (CiiA? + C21A2M, (4) 

with boundary conditions: (/ri(-oo) = 0, \p2{-°°) = 0, A'(-oo) = 1, 1/^1(00) = 1, 1/^2(00) 
1 . a' (00) = 0. The surface energy a„^ is obtained from the substitution of field variables ipi, 1^2, A 
that satisfy the GL eqs (H) into ©. 

It is clear from ([3]) and (U) that the sign of the surface energy is determined by three independent 
dimensionless parameters: ki, K2, k^. If these three parameters are known for a material consid- 
ered, we can then obtain the value of ans/{H'^t'^/Sn) from the substitution of i/rj, iff2, A satisfy 
(HI) into (|3]) and identify the sign of surface energy. If the sign of surface energy is negative, i.e., 
a„.v/(//^;/l/87r) < 0, occurrence of composite vortex is thermodynamically favorable. Otherwise, 
no stable composite vortex can emerge. 
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When the coherence lengths of two components are equal: = ^2, i-C, = 1, it is easily 
verified that equations dH) have the first integral: 
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(5) 

With ([5]), the expression ^ becomes a„s = %^ f^^ dp^^ij/'^ + ^1/^2^ + A' (a' - ijj. The sign 
of surface energy, thus the stability of composite vortex, can be determined by the procedures 
described above. Let us comment on several particular cases: In the case ati » 1, /C2 » 1, the 
first two terms in the integrand can be neglected and the sign of surface energy is always negative 
since A' e [0, 1]. The sign of surface energy is positive in the opposite case ati «; 1, a:2 1. If 
one component is of extreme type-2, while the other component is of extreme type-1, i.e., /ci » 1, 
a:2 «; 1 . The contribution of the first term in the integrand can be neglected and the sign of surface 
energy is positive. 

Let us now study general cases in which there is disparity in coherence lengths between two 
components: ^1 ^2, i-C-j a:^ ^ 1- We then need to solve eqs © with given boundary condi- 
tions numerically . Before a detailed numerical work is undertaken, we first analyze the problem 
qualitatively. We note that in the integrand in ^ only the second term Vd^Pi.iP) contributes the 
negative value to surface energy. The distance scale over which the condensates tends to its expec- 
tation value ¥,0 is of order ~ y(|*F;P) decreases from to a,|¥,oP + bil'Viof/l = -Hl^.^/iSn) 
in the same range. The length scale over which the magnetic field decays is ~ A. The last term of 
integrand in Q ^ (H,., - V x A)^ increases from to ^H-^t in this range. The similar term in the 
surface energy expression when only single superconducting component exists increases from 
to ^HIj(^ in a range ~ Ai. Since A < min(/li, /I2), Hd = (^^(i) + ^ctd)^'^' integral value of the 
last term in integrand in ([2]) is positive and is much larger than the sum of the integral values of 
similar terms in the surface energy expressions when single component exists. We then conclude 
that there is a trend of increase in surface energy for a two-component superconductor comparing 
to the sum of surface energy of the single component cases. The detailed numerical simulations 
below confirm this idea. 

To explore the concrete behavior of the sign of surface energy with three parameters a-i, K2, k^, 
we then search for the numerical solutions of the GL eqs. dH). And we really identified the sign 
change of the surface energy due to the variation of these three parameters. There are three cases: 

Case 1. K\ < II V2, K2 < \l V2, i.e., two components are both of type-1. As we have shown. 
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there is a trend of increase in surface energy for a two-component superconductor comparing to 
the sum of surface energy of the single component cases. Since the sign of surface energy for a 
type-1 material is always positive, we then conclude that the system has positive surface energy 
and there is no thermodynamically stable composite vortex. 

Case 2. > II V2, K2 < 1/ V2, i.e., the first component is of type-2, while the second is 
of type-1. When the third parameter » 1, i.e., ^1 » ^2, thus A\ > ^1/ V2 » ^2 > ^2, 
A = (l//lj -I- l/Aly^^^ X A2 < ^2- Then the integral value of the last term in integrand in ([2]) 
^ (Ret - V X A)^ gains ascendancy over that of the second term y(|*I'i,2p), and the sign of surface 
energy tends to be positive. On the other hand, when /c^ «; 1 , i.e., ^1 «; ^2, penetration depth A 
may fall into the region < A < ^2, and the sign of the surface energy can take negative. We then 
conclude that there is a critical value k^,. at which the surface energy vanishes. As an example, we 
show in Fig. [T]the sign change of surface energy for a two-component superconductor with /ci = 
6.0, K2 = 0.5. It is clear that the sign of surface energy is negative when < k^^ = 0.35. Generally, 
for fixed /Ci > 1 / V2, /C2 < 1 / V2, the critical value k^^ can be determined using numerical method 
as shown above. However, a rough estimate of the upper limit of the critical value k^c can be made 
based on the condition < A < ^2, under which the sign of surface energy can be expected to be 
negative. 



Note that the negative surface energy conduces to the invasion of thermodynamically stable com- 
posite vortices in a two-component material, and these composite vortices assume full responsibil- 
ity for type- 1.5 superconductivity in a type-1 -I- type- 2 superconductor. Then we can regard Eq. ^ 
as a rough criterion for type- 1 .5 superconductivity in the two-component superconductor in which 
type-1 and type-2 condensates coexist. 

Case 3. Ki > 1/ V2, K2 > 1/ V2, i.e., two components are both of type-2. The sign of surface 
energy of the system is then expected to be negative. However, when ki ^ K2~ I, ~ 1, the effect 
of coexistence of the two superconducting components may result in a positive surface energy. In 
Fig. |2]we show the numerical result in the case ki = K2 = 0.9. It is clear to see that there is a sign 
change of the surface energy from negative to positive with increasing parameter k^. This shows 
the the condition conversion for existence and nonexistence of composite vortex. 

At last, we discuss the possible vortex patterns which will occur in a type- 2 -1- type- 2 super- 
conductor after the stability of composite vortex has been confirmed. We concentrate on two 




(6) 
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cases: (i) if the penetration depth is much larger than the coherence lengths of both condensates, 
A » (z = 1,2), two vortices will have their supercurrents overlapping first and repel each other. 
Then, regular vortex patterns like those in an usual type-2 superconductor will occur, (ii) if the 
disparity between the coherence lengths of two components is large, and the penetration depth is 
much smaller than the coherence length of one component, i.e., i «; (or /I «; ^2), the vortex has 
an extended core associated with the condensate *Pi (or ¥2). Interaction between two vortices is 
expected to be short-range repulsive and long-range attractive, similar to that of stable composite 
vortices in the type-l-i-type-2 case. We then expect unusual vortex patterns like those observed in 
MgB2 [Ivl] will occur in these type-2-l-type-2 materials. 

In conclusion, coexistence of two superconducting components in a material may lead to the 
occurrence of stable composite vortex and distinctive response of material to applied magnetic 
field. We have shown that the criterions for stability of composite vortex with normal core in a 
two component superconductor are three independent dimensionless parameters: ki, K2, k^. And 
there can exist stable composite vortex in a type-l-l-type-2 or type- 2 -I- type- 2 superconductor. 

This work was supported by the National Natural Science Foundation of China (No. 
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FIG. 1: Sign change of surface energy for a two component superconductor with K[ = 6.0, K2 = 0.5. It is 
shown that the sign of surface energy changes from positive to negative with decreasing parameter k^. Note 
that the negative surface energy conduces to the invasion of thermodynamically stable vortex. The critical 
value of at which surface energy vanishes is k^c = 0.35. 



8 



0.2 



0.1 - 



00 




FIG. 2: Sign change of surface energy for a two component superconductor with ki - K2 - 0.9. The critical 
value of at which surface energy vanishes is k^c = 0.67. 
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